Precise asymptotics for a variable-range hopping model 



Bernhard Mehlig 1 ) and Michael Wilkinson 2 ) 

1 Department of Physics, Goteborg University, 41296 Gothenburg, Sweden 
2 Faculty of Mathematics and Computing, The Open University, Walton Hall, 
Milton Keynes, MK7 6AA, England 

For a system of localised electron states the DC conductivity vanishes at zero temper- 
ature, but localised electrons can conduct at finite temperature. Mott gave a theory for 
the low-temperature conductivity in terms of a variable-range hopping model, which is hard 
to analyse. Here we give precise asymptotic results for a modified variable-range hopping 
model, proposed by S. Alexander [Phys. Rev. B, 26, 2956 (1982)]. 



§1. Introduction 

If a system has localised electron states the DC conductivity must be zero at 
zero temperature, but localised electrons can conduct at finite temperatures. Mott 1 ) 
proposed that the low-temperature dependence of the DC conductivity is 

(1-1) 

in d spatial dimensions (variable-range hopping law). 

The variable-range hopping process is described by specifying transition rates 
for an electron to hop from one localised site to another. The standard model 
is illustrated in one spatial dimension in figure [TJ It is of the form of a resistor 
network, with nodes labelled by n = 1, . . . , N. The conductances T nm are given by 
the transition rates 

r nm = r exp (^ J E n\ + \Ern\+^E n -E m \ ^ \ n _ m | ^ 

where E m are random on-site energies (measured relative to the Fermi energy). 
Here the exponential factor accounts for the temperature dependence of the rate 
constants and the factor el n_m l models the decrease of the matrix elements as a 
function of distance between the sites. The matrix elements are assumed to decrease 
exponentially as a function of hopping distance, with parameter e = exp(— a) (where 
a is the inverse localisation length). In the low-temperature limit the rate constants 
differ by orders of magnitude, so that the conduction of the chain is determined by 
a single path, and the conductivity is the harmonic mean of the dominant transition 
rates. Transitions have a higher rate if they require a small increase in the energy 
of the electron, which favours finding a conduction path involving low-energy states. 
There is thus a competition between hopping to a nearby site, which is favoured by 
large transition matrix elements, and finding an energetically favoured site, which 
may be some distance away. 
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Fig. 1. One-dimensional resistor-network model for a one-dimensional variable-range hopping prob- 
lem. The conductances _T nm are given by (11-2[1 . E n are random on-site energies. 



Despite the appealing simplicity of this physical picture, there are no precise 
asymptotic expressions and few significant results on variable-range hopping. Am- 
begaokar et al?' have employed the fact that the transition rates between localised 
electron states in disordered solids fluctuate wildly to map the variable-range hop- 
ping problem onto a percolation problem. This explains the exponent in (jl-lj) for 
d > 1 and provides an estimate of To. Kurkijarvi 3 ) showed that the DC conduc- 
tivity in infinitely long disordered chains is simple activated conduction, following 
the Arrhenius law: o~(T) ~ exp(— Tq/T). This result was later refined by Raikh and 
Ruzin. 4 ) In one dimension, Mott's law (|l-ip is expected to be valid for sufficiently 
short chains only. 

Alexander 5 ) introduced a simplified model for variable-range hopping. In this 
paper we discuss the physical significance of his model and analyse it using an al- 
ternative method. We obtain a precise asymptotic expression for the conductivity 
according to this model. 

§2. The Alexander model 

Finding the true conduction path is a difficult problem. Here, following Ref. [5], 
we adopt a simpler approach. We construct a conduction path, visiting a sequence 
of sites by the following scheme. Given the ith site in the sequence, at position n 
in the lattice, we take the site i + 1 to be the site with position m which has the 
largest transition rate T nm [given by equation (|l-2p ] , considering only sites lying to 
the right. (That is, we maximise T nm with respect to m, subject to m > n). This 
approach will not usually produce the optimal conduction path, for which a non- 
optimal transition at one step may be more than compensated by a lower overall 
resistance. However, the model has the appealing feature that its asymptotics can 
be determined precisely, as we now show. 

We can simplify the model further by assuming that the rate is F nm = ^"-H^^-V* 
where uj n is a random number and t a small parameter. In the following we assume 
that the uo n are independently, identically uniformly distributed in [0, 1]. 
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§3. Method of solution 

Because the values of r nm differ by orders of magnitude we can approximate the 
largest rate r from a given site to some other site by the sum over all rates 

oo 



r = Y j e n x n (3-1) 

n=0 

where x n = e~ U)n l t . The probability density for the random variables x n is 

p t x \ = { t / x for exp(-t" 1 ) < x < 1 
y otherwise. 

The conductivity of the system is obtained from the diffusion constant D, which is 
given by D = —JAN/Ap, where Ap is the change in occupation probability over 
AN sites. We have Ap = JJ^i^l ~ 1 = J(r~ 1 )AN stcp , where AN step is the number 
of maximum-rate hops in the process described above. If the maximum-rate hops 
have average length (An), then AN step = AN /(An), and the diffusion constant is 

D = (r- 1 )- 1 (An) . (3-3) 

The quantity (An) may be determined by considering 

^aiogT _ Y, n ne n x n 



(34) 



This will almost always be approximately equal to the value of n for which the rate 
e n x n is maximal. Thus we may write 

(An) =e^ (log r) (3-5) 

and we have the ingredients necessary to determine D from the probability density 
pr of r. 

§4. Distribution of T 

The distribution of r may be obtained from p(x) [the distribution of the vari- 
ables x n in (|3Tp ]. by application of the convolution theorem. Because J 1 is a sum of 
independent random variables, its probability density is the convolution of the prob- 
ability density of each component, namely p(x/e n ). Correspondingly, the generalised 
Fourier transform pr(k) of the required probability density is a product: 

oo 

pr(k) = l[p(ke n ) (4-1) 

71=0 

where the generalised Fourier transform is defined by 

/•oo 

p(k) = dx exp(ikx)p(x) . (4-2) 
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Fig. 2. (Left) shows results of numerical simulations for diffusion constant D as a function of t 
(o) for e = 0.2, compared to the exact asymptotic theory, equation (|7-2[) . (Right) shows S(a), 
equation (|7-3[) . as a function of a (solid line). Also shown is the approximate result S(a) ~ —a/2 
from Ref. [5] (dashed line). The arrow indicates the nearest-neighbour limit S = — 1. 



Setting u = kx we obtain 



/A; 
dn n -1 exp(ru) 

k exp(— 1/t) 

= *[J5i(-ifcexp(-l/<) - Ei(-ifc)] 
where E\[v) is the exponential integral 



(4-3) 



cxp( 



(4-4) 



From the asymptotic form for E\(v) we find p(0) = 1 (which is required by normal- 
isation). Further we find \p(k)\ < t for k > exp(l/t), and that for 1 <C k <C exp(l/i) 
we can apply the following approximation: 



~ 1 — t log 



jt + ii7Tt . 



(4-5) 



We now determine pr{k). In the case where \k\ <C 1 all of the factors are close to 
unity such that pr(k) ~ 1 in that region. In the case where \k\ 3> exp(l/i), the 
product is negligible, and we approximate p~r(k) ~ in that region. Finally, when 
1 <C k -C exp(l/t), most of the factors in the infinite product are very close to unity, 
and we can replace the infinite sum by a finite number N* of terms. This number is 
determined by the value TV* of the index such that the argument of the corresponding 
term, e N * k is equal to unity. Thus N* is the integer part log fe/log(l/e) = logk/a 
and we have 

Int[(logfc)/a] 

Pr( k ) ~ ]1 {l-tloge n k-~ft+±Trit) . (4-6) 

n=0 
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Evaluating the product we obtain an approximate expression for pr(k), valid for 
t <C 1 and 1 < k < exp(l/i): 



r 



1 + at - 7t + i7rt/2 



pr(fc) ~ (at) 1+Q " llogfe ,) nTTT^r , r-r4 . (4-7) 



r 



/l -tlog/c-7t + ivrt/2\ ' 
V at / 



We need to determine the asymptotic form for small values of t. Note that x = log k 
varies between and t" 1 so that xt is of order unity. Using 

r(z + a) ~ V2^ e - z z~^ +z+a (4-8) 

for large z and taking z = (1 — xt)/(at), we obtain 

x , 1 1 — xt 7 ivr N , , 
+ ( + ) log(l - xt) 

(4-9) 

The calculation of D requires two averages: (r~ l ) and (logT). Averages can be 
calculated in terms of the generalised Fourier transform of the probability density: 

/CO -1 />oo 

dr P (r)/(r) = — J dk ~ Pr (k)~f{k) (4-io) 

where f(k) is the generalised Fourier transform of the weight function f(r). 

§5. Harmonic average of r 

We have 

roo 

(r- 1 ) = / dk lm[p r (k)] (5-1) 
J o 

since the generalised Fourier transform of 1/.T is 7risgn(fc). We evaluate the integral 
in the saddle-point approximation. To this end let y = xt = t log k. Note that the 
contribution to integrals with respect to k from the interval [0, 1] is negligible and is 
not considered. We thus have 

1 f 1 

(J- 1 )--/ dye^ ls ^e^ + ^) 1 °g( 1 -f) S m^(y) (5-2) 
t Jo 



with 



S{y)=y-y-- l —^-\og{l-y), (5-3) 
a a 

Hy) = TT-iogU - y) ■ 

la 



The saddle-point equations are: 
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d 2 S* 



c 



It follows that 



dy 2 a 

-l + a + e"" 



(r- 1 ) ~y^ e - a -Te at . (5-5) 

§6. Average of log r 



To calculate (log-T), we make use of the relation 



(lo g r}~^- [\y( y + 1 t)e {at) 1 - R ^e ( H^) 1 °s(i-y) sin 0( y ). (g . 2) 



(logr) = - f°°dk 1 ° gl ^ l+7 Im[ j 5 r (A;)] . (6-1) 
n Jo k 

From equation (|6-ip we find 
1 2 r 

with 

R(y) = -y - (I - y)\og{\ - y) (6-3) 
and <j) as defined in (|5-3p . The integral is approximately 



(logr)--i-^- [\yy(y + 1 t)e^ lR ^ . (64) 



t 2 vr 2a _ 

Approximating i?(y) = —y 2 /2, we obtain 

(logr)~-(7T/2) 1 /2( a / t )l/2. ( 6 . 5 ) 

§7. Results 

Using (|3-5p it follows that 



Taking equations (|3-3|) . (j5-5jl . and (|7-ip together we finally obtain: 

J) ^ J_ e a+7 e [l-a~exp(-a)]/(at) ^.2) 

4a 

This is the main result of this paper. It describes the exact asymptotics of the 
diffusion constant in Alexander's model, exhibiting Arrhenius behaviour, with action 
(excitation energy) given by 

S(a) = 1 " a ~ eXp( " a) (7-3) 
a 
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which is plotted in Fig. [2j The limit 5(0) = —a/2 (dashed line in Fig. [2]) was 
previously obtained in Ref. [5]. In the opposite limit (a — ► oo) the action converges 
to —1, describing activated transport in a conventional nearest-neighbour model. 

In addition to the precise form (|7-3j) of the action, our result (|7-2|) also gives the 
pre-exponential prefactor exactly. We know of no variable-range hopping model for 
which the pre-exponential factor has been obtained precisely. The precise asymptotic 
expression for D allows for comparison with computer simulations. The simulations 
were performed as follows: we sampled the rate r for the optimal hop from a given 
point to the right, as well as the corresponding length An. The diffusion constant 
was then calculated from (|3-3|) . In figure [2j equation (|7-2p is shown together with 
results of computer simulations for e = 0.2 (corresponding to a ~ 1.62). We observe 
excellent agreement between the simulations and the precise asymptotic theory. 
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